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Bulk  conduction  and  relaxation  of  the  [(Zr02)i-x(Ce02)x]o.92(Y203)o.o8  (0<x<l)  solid  solutions  were 
studied  using  impedance  spectroscopy  at  intermediate  temperatures  (200-500  °C).  The  bulk  conductivity 
as  a  function  of  x  shows  a  “V-shape”  variation  which  is  a  competitive  effect  of  the  defect  associates  and 
the  lattice  parameter.  In  the  Zr02-rich  region  (x<  0.5)  Ce02  doping  increases  the  concentration  of  defect 
associates  which  limits  the  mobility  of  the  oxygen  vacancies;  in  the  Ce02-rich  region  (x  >  0.5)  the  increase 
of  x  increases  the  lattice  parameter  which  enlarges  the  free  channel  for  oxygen  vacancy  migration.  Further 
analysis  indicates  the  ionic  radius  of  the  tetravalent  dopant  determines  the  composition  dependence  of 
the  ionic  conductivity  of  the  solid  solutions.  When  doping  YSZ  with  other  tetravalent  dopant  with  similar 
ionic  radius  with  Zr4+,  e.g.,  Hf4+,  such  “V-shape”  composition  dependence  of  the  bulk  conductivity  cannot 
be  observed. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Oxide  ion  conductors  with  cubic  fluorite  structure  are  impor¬ 
tant  electrolyte  materials  for  solid  oxide  fuel  cells  (SOFCs)  [1,2], 
oxygen  sensors  [3]  and  oxygen  pumps  [4].  High  ionic  conductivity 
is  achieved  in  such  a  structure  when  the  host  cations  are  replaced 
by  lower  valence  cations,  which  generates  oxygen  ion  vacancies  to 
compensate  the  charge  missing  and  to  act  as  charge  carriers  [5]. 
The  oxygen  ions  overcome  an  energy  barrier  to  hop  to  neighbour¬ 
ing  vacant  sites  and  move  diffusively  to  fulfil  long-range  transport 
in  the  material  [6]. 

Among  numerous  oxide  ion  conductors,  Zr02-  and  Ce02-based 
ceramic  attract  the  most  interest  and  have  been  extensively  studied 
[7-15].  They  possess  many  unique  physical  properties  which  make 
them  promising  materials  for  practical  or  potential  application  in 
technologically  important  devices.  For  example,  yttria-stabilized 
zirconia  (YSZ)  is  the  conventionally  employed  electrolyte  material 
in  SOFCs  because  of  its  good  mechanical,  chemical  and  electrolytic 
properties  [16].  Doped-ceria  has  higher  ionic  conductivity,  good 
thermodynamic  stability  [17]  and  good  catalytic  activity  [18], 
therefore  it  is  a  prime  candidate  for  the  intermediate  temperature 
SOFCs. 

The  Zr02-Ce02-Y203  ternary  solid  solution  which  keeps  the 
cubic  fluorite  structure  is  also  a  material  of  interest  both  from  the 
practical  side  and  the  fundamental  side.  From  the  practical  side, 
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the  idea  of  a  double-layer  electrolyte  with  YSZ  on  the  anode  side 
and  doped  Ce02  on  the  cathode  side  has  been  proposed  for  SOFC 
[19,20].  The  double-layer  electrolyte  combines  the  advantages  of 
YSZ  and  Ce02:  it  suppresses  the  electronic  conduction  caused  by 
the  reduction  of  Ce4+  at  the  anode  side  and  retains  the  high  ionic 
conductivity  of  doped-Ce02  [21].  However,  the  reaction  of  Zr02 
and  Ce02  leads  to  the  formation  of  the  Zr02-Ce02-Y203  ternary 
solid  solution  between  the  two  layers.  Therefore,  an  investigation 
of  the  ionic  conductivity  of  the  ternary  solid  solution  is  necessary. 
On  the  fundamental  side,  it  is  a  system  to  study  the  homovalent 
doping  effect  on  the  ionic  conductivity  of  the  solid  solutions  if  the 
concentration  of  Y203  is  constant. 

Several  studies  on  the  electrical  properties  of  the 
Zr02-Ce02-Y203  solid  solutions  have  been  reported.  How¬ 
ever,  most  of  the  studies  are  focused  on  the  mixed  conduction 
of  the  solid  solutions  caused  by  the  reduction  of  Ce4+  to  Ce3+ 
at  various  oxygen  partial  pressures  and  at  high  temperatures 
(>800 °C).  For  example,  Cales  and  Baumard  [22]  studied  the  total 
conductivity  of  the  ternary  solid  solutions  in  a  large  range  of 
oxygen  partial  pressures  from  1000  to  1400  °C.  They  discussed 
the  mixed  conduction  and  the  defect  structure  caused  by  the 
reduction  of  Ce4+  to  Ce3+,  and  determined  the  concentrations 
of  the  various  charge  carriers  by  thermogravimetry  and  the 
measurement  of  the  magnetic  susceptibility.  Arashi  et  al.  [23] 
studied  the  electrical  conduction  of  the  Zr02-Ce02-Y203  solid 
solutions  at  even  higher  temperatures  (>1200  °C).  They  employed 
the  electron  blocking  method  to  separate  the  ionic  and  electronic 
conductivities  and  found  that  the  reduction  of  Ce4+  to  Ce3+  not 
only  influences  the  relationship  between  ionic  conductivity  and 
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the  oxygen  partial  pressure,  but  also  influences  the  electronic 
contribution  to  the  total  conductivity.  Ananthapadmanabhan  et  al. 
[24]  studied  the  composition  dependence  of  the  electrical  conduc¬ 
tivity  in  [(Zr02)i_x(Ce02)x]0.9(Y203)o.i  solid  solutions  from  ~700 
to  1300  °C  in  air  and  they  proposed  a  simple  model  for  the  oxygen 
path  length  to  explain  the  observed  composition  dependence  of 
the  conductivity. 

The  reduction  of  Ce4+  at  high  temperatures  makes  the 
Zr02-Ce02-Y203  solid  solution  a  complicated  system  to  study  the 
effect  of  tetravalent  doping  on  the  ionic  conductivity.  On  the  other 
hand,  most  of  the  current  studies  on  the  electrical  properties  of  the 
Zr02-Ce02-Y203  solid  solutions  were  studied  by  DC  four-probe 
measurement.  The  conductivity  obtained  from  the  dc  measurement 
is  an  overall  conductivity  from  both  grain  (bulk)  and  grain  bound¬ 
ary.  In  polycrystalline  ceramic  materials,  the  grain  boundary  has  a 
large  influence  on  the  electrical  conductivity  because  of  impurity 
segregations  or  the  space  charge  effect  [25].  For  studying  the  intrin¬ 
sic  property  of  the  ternary  solid  solution,  the  contribution  of  grain 
boundary  to  the  overall  conductivity  should  be  excluded. 

Therefore,  in  this  paper,  the  electrical  properties  of  the 
[(Zr02)i_x(Ce02)x]o.92(Y203)o.o8  (0<x<l)  were  studied  in  an 
intermediate  temperature  range  (200-500  °C)  by  impedance  spec¬ 
troscopy.  The  relative  low  temperatures  were  chosen  to  avoid  the 
reduction  of  Ce4+  and  to  ensure  the  electrical  conduction  in  the  solid 
solutions  is  purely  ionic.  Impedance  spectroscopy  was  employed  to 
separate  the  bulk  and  the  grain  boundary  contributions  to  the  elec¬ 
trical  conductivity.  The  purpose  of  this  work  is  to  study  the  effect  of 
homovalent  doping  on  the  ionic  conductivity  of  the  solid  solutions. 


2.  Experiment 
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Fig.  1.  X-ray  diffraction  patterns  of  the  solid  solutions.  The  inset  is  the  lattice  param¬ 
eter  as  a  function  of  the  mole  ratio  of  Ce02  (x)  in  [ (Zr02)i_x(Ce02)x  ]o.92(Y2 03)0.08- 


were  normalized  by  a  geometric  factor  oi  =  AtlnD 2,  where  t  and  D 
denote  to  sample  thickness  and  diameter,  respectively. 

3.  Results 

3.1.  Phase  composition  and  lattice  parameter 

Fig.  1  shows  the  X-ray  diffraction  patterns  of  the  solid  solutions, 
which  confirms  that  all  the  samples  are  single  phase  with  cubic 
structure.  The  lattice  parameters  of  the  samples  are  shown  in  the 
inset  of  Fig.  1,  where  linear  relationship  between  the  lattice  param¬ 
eter  and  the  mole  ratio  of  Ce02  can  be  observed,  indicating  it  follows 
the  Vegard’s  rule  for  solid  solutions. 


The  samples,  formulated  as  [(Zr02)i_x(Ce02)x]o.92(Y203)o.o8 
(0  <  x  <  1 ),  were  synthesized  by  a  solid  state  reaction  method.  Zr02 
(99.6%,  PI-KEM,  UK),  Y203  (99.99%,  PI-KEM,  UK)  and  Ce02  (99.95%, 
PI-KEM,  UK)  powders  were  used  as  the  starting  materials.  Appropri¬ 
ate  amount  of  each  powders  were  weighed,  mixed  and  ball-milled 
in  2-propanol  for  24  h  using  zirconia  balls  as  grinding  media. 
The  resulting  mixtures  were  subsequently  dried  in  air  overnight, 
and  then  ground  by  mortar  and  pestle  and  passed  through  a  45 
micro  sieve.  The  final  mixtures  were  cold  pressed  into  tablets 
under  a  uniaxial  pressure  of  100  MPa  and  then  sintered  at  1500°C 
for  10  h  in  air.  The  other  series  of  solid  solutions,  formulated  as 
[(Zr02)i_x(Hf02)x]o.87(Y203)o.i3  (0<x<0.8),  were  synthesized  by 
the  same  powder  mixing  method,  and  sintered  at  1600  °C  for  4  h. 

Phase  compositions  of  the  solid  solutions  were  identified  by 
X-ray  diffraction  (XRD,  Philips  X’Pert)  method  using  Cu  Ka  radia¬ 
tion.  The  measurements  were  performed  on  the  sample  surfaces 
with  a  step  scanning  mode  (step  size  of  0.05°)  at  a  rage  of 
0.1  min-1.  The  densities  of  the  sintered  specimens  were  measured 
by  Archimedes’  method.  The  porosities  of  all  the  samples  are  within 
10%.  Microstructures  of  the  samples  were  observed  by  scanning 
electron  microscope  (SEM,  Philips  XL30).  The  grain  sizes  of  the  sam¬ 
ples  were  estimated  from  the  SEM  images,  and  they  were  found  to 
be  around  10  pan  for  all  the  compositions. 

Electrical  properties  of  the  solid  solutions  were  obtained  from 
ac  impedance  spectroscopy  measurements  using  a  Solatron  SI  1 255 
HF  frequency  response  analyser  coupled  with  a  Solatron  1296 
Delectric  Interface  (Solartron,  UK).  Silver  paint  was  coated  on  the 
polished  surfaces  of  the  samples  and  fired  at  500  °C  for  1  h  in  air  to 
serve  as  electrodes.  During  impedance  measurements,  an  ac  volt¬ 
age  of  0.1  V  was  applied  to  the  sample  over  a  frequency  range  from 
0.1  to  107Hz  at  various  temperatures.  Equivalent  circuit  fittings 
of  the  measured  impedance  spectra  were  carried  out  using  Zview 
Impedance  Analysis  software  (Scribner  Associates,  Inc.,  Southern 
Pines,  NC).  For  comparison  convenience,  all  the  impedance  spectra 


3.2.  Typical  impedance  spectra  of  the  solid  solutions 

First  of  all,  one  composition  of  the  solid  solutions 
[(Zr02)0.6(Ce02)oi4]o  .92(Y2O3)0.08  (x  =  0.4)  was  selected  to  present 
the  typical  ac  impedance  behaviours  of  the  solid  solutions.  Fig.  2 
shows  typical  impedance  spectra  (Nyquist  plots)  of  the  solid 
solutions  at  various  temperatures.  Similar  to  the  well-known 
impedance  spectroscopy  of  YSZ,  two  semicircles  are  displayed 
on  the  Nyquist  plot  of  the  solid  solution,  from  left  to  right  (high 
frequency  to  low  frequency),  representing  the  response  from  grain 
(bulk)  and  grain  boundary,  respectively.  The  measured  impedance 
spectra  can  be  fitted  by  an  equivalent  circuit  of  two  parallel-aligned 
resistance-constant  phase  element  (R-CPE)  in  series  connection,  as 


p'  (ftm) 

Fig.  2.  Typical  Nyquist  plots  of  the  solid  solutions  at  various  temperatures  (x  =  0.4 
in  this  case).  The  red  open  circles  are  experimental  values,  while  the  black  solid 
lines  are  the  equivalent  circuit  fitting  results.  The  inset  figure  on  the  top-left  is  the 
equivalent  circuit,  and  the  top-right  inset  figure  shows  the  Arrhenius  plots  of  the 
bulk  and  grain  boundary  conductivity. 


F.  Yang  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  4943-4949 


4945 


Frequency  (Hz) 


0.020 
0.016 
^  0.012 
S  0.008 
0.004 
0.000 

101  102  1  03  1  04  1  05  1  06 

Frequency  (Hz) 


Fig.  3.  Bode  plots  (imaginary  impedance  vs.  frequency)  of  the  solid  solution  (x  =  0.4) 
at  various  temperatures.  The  top-left  inset  figure  shows  the  Arrhenius  plot  of  the 
bulk  relaxation  frequency.  The  top-right  inset  figure  shows  the  scaling  behaviour  of 
the  imaginary  impedance. 


Fig.  4.  Frequency  dependence  of  the  imaginary  electric  modulus  (M")  of  the  solid 
solution  (x  =  0.4)  at  various  temperatures.  The  top-left  inset  figure  shows  the  Arrhe¬ 
nius  plot  of  the  bulk  relaxation  frequency.  The  top-right  inset  figure  shows  the 
scaling  behaviour  of  the  imaginary  electric  modulus. 


illustrated  by  the  top-left  inset  in  Fig.  2.  The  bulk  and  grain  bound¬ 
ary  resistances  were  respectively  obtained  from  the  equivalent 
circuit  fitting  results  and  subsequently  converted  to  conductivity 
by 


° bulk  ~ 


1  t 
Rbulk  A 


and 


°£  = 


2  Cbulk  2 
rgb  Cgb  A 


(1) 


(2) 


The  grain  boundary  conductivity  is  2  orders  of  magnitude  lower 
than  the  bulk  conductivity,  which  is  caused  either  by  the  impurity 
segregation  at  the  grain  boundary  or  the  space  charge  effect  due  to 
the  oxygen  vacancy  depletion  near  the  grain  boundary  region  [25]. 
The  temperature  dependence  of  the  bulk  and  the  grain  boundary 
conductivity  is  shown  in  the  top-right  inset  figure  in  Fig.  2.  Both  the 
bulk  and  the  grain  boundary  conductivity  shows  a  typical  thermally 
activated  process  that  it  increases  with  an  increase  in  the  temper¬ 
ature  and  its  dependence  on  the  temperature  obeys  the  Arrhenius 
law  in  the  form  of 


crif  o  exp 


(3) 


where  i  represents  bulk  or  grain  boundary,  cr0  is  a  pre-exponential 
factor,  1<b  is  the  Boltzmann  constant,  T  is  the  absolute  tempera¬ 
ture  and  Ea  is  the  activation  energy.  The  activation  energies  can 
be  obtained  from  the  slope  of  the  Arrhenius  plots,  with  values 
of  1.08  ±0.01  eV  for  bulk  conduction  and  1.09  ±  0.05  eV  for  grain 
boundary  conduction. 

Fig.  3  shows  the  Bode  plots  (imaginary  impedance  vs.  frequency) 
of  the  selected  composition  (x  =  0.4)  at  various  temperatures.  The 
two  peaks  on  the  Bode  plot  correspond  to  the  two  semicircles  on 
the  Nyquist  plot  and  the  peak  height  is  proportional  to  the  resis¬ 
tance  of  each  response  (diameter  of  each  semicircle).  Because  the 
grain  boundary  resistance  is  much  smaller  than  the  bulk  resistance 
in  this  sample  (see  Fig.  2),  the  low  frequency  peak  corresponding 
to  the  grain  boundary  response  is  not  obviously  distinguished  on 
the  Bode  plot,  but  this  does  not  influence  the  discussion  since  we 
mainly  focus  on  the  high  frequency  response  (bulk  conduction).  The 
frequency  at  which  the  peak  reaches  its  maximum  is  defined  as  a 
relaxation  frequency  (fZ").  As  shown  in  the  top-left  inset  in  Fig.  3, 


the  temperature-dependence  of fz»  follows  the  Arrhenius  law  given 
by 

/z»=/oexP  (-gf),  (4) 

where /o  is  the  pre-exponential  factor,  EZa  is  the  activation  energy. 
The  activation  energy  obtained  from  fz»  is  1 .05  ±  0.02  eV.  The 
top-right  inset  in  Fig.  3  shows  the  scaling  results  at  different  tem¬ 
peratures  of  the  imaginary  impedance  where  p^ax  and  fz»  max  are 
used  as  the  scaling  parameters  for  p"  and/,  respectively.  It  can  be 
seen  all  the  imaginary  impedance  collapses  into  one  master  curve, 
indicating  the  conduction  mechanism  remains  unchanged  within 
this  temperature  range  [26]. 

The  frequency  dependence  of  the  imaginary  electric  modulus 
(M")  at  various  temperatures  is  shown  in  Fig.  4.  In  the  modulus  plot, 
only  one  peak  corresponding  to  the  bulk  is  displayed  because  of  its 
much  smaller  capacitance  value  than  that  of  grain  boundary.  The 
peak  maximum  M(^ax  slightly  increases  with  increasing  tempera¬ 
ture  indicating  a  weak  temperature  dependence  of  the  dielectric 
constant  of  the  sample.  Similarly,  the  relaxation  frequency  at 
which  the  imaginary  electric  modulus  reaches  its  maximum  value, 
defined  as  /m",  shifts  to  higher  frequency  range  with  increas¬ 
ing  temperature,  and  its  temperature-dependence  also  obeys  the 
Arrhenius  law,  as  shown  in  the  top-left  inset  of  Fig.  4.  The  activa¬ 
tion  energy  (EMa)  value  is  1.07  ±  0.03  eV  which  is  similar  to  EZa.  In 
a  scaled  coordinate,  as  shown  in  the  top-right  inset  in  Fig.  4,  all  the 
curves  collapse  into  a  single  master  curve,  indicating  the  relaxation 
describes  the  same  mechanism  at  various  temperatures  [27]. 

In  Fig.  5,  the  imaginary  impedance  (p")  and  the  imaginary  mod¬ 
ulus  (M")  were  plotted  as  a  function  of  frequency.  The  peaks  are  not 
overlapped:  the  peak  position  for  M"  shifts  to  a  higher  frequency 
region  compared  to  the  p"  peak.  According  to  Ref.  [28],  the  overlap¬ 
ping  peak  position  of  imaginary  impedance  and  imaginary  modulus 
is  evidence  of  delocalized  or  long-range  relaxation.  Therefore,  the 
slight  separation  of  the  two  peaks  (the  peak  position  of  M"  is  around 
half  magnitude  higher  than  that  of  p”)  suggests  the  components 
from  both  long-range  and  localized  relaxation  [29]. 

3.3.  Composition  dependence  of  the  bulk  conductivity 

Fig.  6  shows  the  composition  dependence  of  the  bulk  conductiv¬ 
ity  of  the  solid  solutions  at  various  temperatures.  With  the  increase 
in  the  Ce02  mole  ratio,  the  conductivity  first  decreases,  reaching  a 
minimum  around  x  =  0.5,  and  increases  afterward.  The  “V-shape” 
variation  of  the  bulk  conductivity  as  a  function  of  the  mole  ratio  of 


4946 


F.  Yang  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  4943-4949 


Fig.  5.  Frequency  dependence  of  the  imaginary  impedance  (p")  and  the  imaginary 
electric  modulus  (M")  of  the  solid  solution  (x  =  0.4)  at  300  °C. 


Ce02  indicates  that  the  ternary  solid  solution  deteriorates  the  elec¬ 
trical  conductivity  of  (ZrO2)0.92(Y2O3)0.08  and  (CeO2)0.92(Y2O3)0.08 
binary  systems.  The  inset  in  Fig.  6  shows  the  relaxation  frequen¬ 
cies  obtained  from  imaginary  impedance  and  imaginary  modulus 
and  their  dependence  on  the  composition  of  the  solid  solutions  at  a 
chosen  temperature.  It  is  noticeable  that  fz«  is  lower  than  fM"  for  all 
the  compositions,  indicating  the  existence  of  both  long-range  and 
localized  relaxations  in  the  solid  solutions.  On  the  other  hand,  the 
composition  dependence  of  fZ"  and  fM"  has  a  similar  trend  as  the 
bulk  conductivity. 

Fig.  7  shows  the  temperature  dependence  of  the  activation 
energy  for  the  bulk  conduction.  The  activation  energy  has  a  decreas¬ 
ing  trend,  although  not  linearly,  with  an  increase  in  x.  Slight 
decrease  of  the  activation  energy  is  observed  when  x<0.5,  while 
a  rapid  drop  of  the  activation  energy  occurs  when  x  >  0.5. 


4.  Discussions 

4.1.  Relationship  between  the  relaxation  frequencies  and  the 
conductivity 

On  the  Bode  plot  (imaginary  impedance  vs.  frequency),  a 
relaxation  frequency  is  obtained  where  the  imaginary  impedance 
reaches  its  maximum  value.  In  an  ideal  case  of  an  equivalent  circuit 
of  a  resistor  and  a  capacitor  in  parallel  connection,  the  imaginary 


Fig.  6.  Composition  dependence  of  the  bulk  conductivity  of  the 
[(Zr02)i-x(Ce02)x]o.92(Y2 03)0.08  (0<x<l)  solid  solutions  at  various  tempera¬ 
tures.  The  inset  figure  shows  the  composition  dependence  of  relaxation  frequencies 
obtained  from  imaginary  modulus  (fMv )  and  imaginary  impedance  (fzn )  at  300  °C. 


Fig.  7.  Composition  dependence  of  the  activation  energy  for  the  bulk  conduction. 
The  hollow  triangles  are  the  experimental  values  obtained  from  the  Arrhenius  plot 
of  the  bulk  conductivity.  The  red  dot  is  the  trend  line.  (For  interpretation  of  the 
references  to  color  in  this  figure  caption,  the  reader  is  referred  to  the  web  version 
of  the  article.) 


impedance  is  expressed  as: 
coR2C 


- Z "  = 


1  +  co2R2C2' 


(5) 


where  co  is  the  angular  frequency,  R  is  the  resistance  and  C  is  the 
capacitance.  Maximum  of  -Z"  is  obtained  when  co  =  1  IRC.  The  bulk 
resistance  and  capacitance  are  determined  by  its  conductivity  (cr), 
dielectric  constant  (e)  and  its  geometries  (t  and  A  as  defined  in  Eq. 
(2))  by 


R=  M 
g  A 


(6) 


and 


Therefore  the  relaxation  frequency  fz"  can  be  expressed  by 


fz"  = 


1  a 
2n  s ' 


(7) 

(8) 


In  Zr02  based  ceramics,  the  dielectric  constant  is  usually  insensitive 
to  the  composition  of  the  material  [30,31]  As  a  result,  fz»  is  only 
determined  by  the  conductivity  (Tm"  has  the  same  expression  as 
fz"  in  Eq.  (8)).  A  direct  proportional  relationship  can  be  established 
between  the  relaxation  frequency  and  the  conductivity. 

In  the  real  case,  the  capacitance  response  is  not  a  pure  capacitor, 
thus  the  expression  for  fzn  will  be  different  from  Eq.  (8).  But  for 
bulk  impedance  (high  frequency  response  on  the  Nyquist  plot),  the 
deviation  from  a  capacitor  is  small.  Therefore,  Eq.  (8)  can  still  be 
used  to  describe  the  relationship  between  fz»  and  a.  Comparing 
the  relaxation  frequency  can  be  a  simple  method  to  rank  the  bulk 
conductivity  of  the  solid  solutions. 


4.2.  Origin  of  the  localized  relaxation/ conduction 

When  Zr02  or  Ce02  is  doped  with  Y203,  oxygen  vacancies  are 
generated  via  the  defect  reaction: 

Y203^2Y'm+Vo  +  30x0, 

where  M  represents  Zr  or  Ce.  However,  due  to  the  coulomb  inter¬ 
action,  some  of  the  oxygen  vacancies  may  bind  to  the  dopant  ions 
and  form  the  defect  associates  (Y'M  -  Vq  ).  It  is  usually  accepted  the 
orientation  of  the  defect  associates  under  an  applied  electric  field 
contributes  to  the  localized  relaxation  [32].  On  the  other  hand,  the 


F.  Yang  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  4943-4949 


4947 


•  9  9  •  • 

9  9  9 

9  9  9  9  9 

9  9  9  9 

9  9  9  9  9 

9  9  9  9 

9  9  9  9  9 

9  9  9  9 

9  9  9  9  9 

Fig.  8.  A  schematic  top  view  of  the  possible  sites  of  oxygen  vacancies  around  a 
dopant  ion.  Numbers  1,  2  and  3  represent  INN,  2NN  and  3NN  oxygen  ions  with 
respect  to  the  dopant  (Y)  ion. 

free  oxygen  vacancies  also  contribute  to  a  localized  conduction,  as 
discussed  below. 

The  ionic  radii  of  Zr4+  and  Y3+  are  0.084  and  0.1 01 9  nm,  respec¬ 
tively  [33],  while  the  ionic  radius  of  Ce4+  has  different  reference 
values,  varies  from  0.097  [34],  0.101  [35]  to  O.llOnm  [36].  How¬ 
ever,  it  is  clear  that  the  ionic  radius  of  Y3+  is  close  to  Ce4+  but  it  is 
significantly  larger  than  Zr4+.  Density  functional  theoretical  studies 
by  Andersson  et  al.  [37]  reported  that  the  formation  of  an  anionic 
vacancy  is  correlated  to  the  ionic  radius  of  the  dopant  ion.  With 
respect  to  the  dopant  ion,  there  are  three  positions,  named  the 
first  nearest  neighbour  ( 1 NN),  second  nearest  neighbour  (2NN)  and 
third  nearest  neighbour  (3NN),  for  the  oxygen  vacancy  to  locate, 
as  illustrated  in  Fig.  8.  INN  position  is  favourable  for  the  oxygen 
vacancy  when  the  dopant  ion  has  an  ionic  radius  closer  to  the  host 
cation,  and  the  2NN  position  is  favourable  when  the  dopant  ion 
has  an  ionic  radius  considerably  larger  than  the  host  cation.  There¬ 
fore,  when  Y3+  ions  substitutes  Zr4+  ions,  oxygen  vacancy  locates 
at  the  2NN  positions  (close  to  Zr  centre);  when  Y3+  ions  substitutes 
Ce4+  ions,  oxygen  vacancy  locates  at  the  INN  positions  (close  to  Y 
centre). 

Dholabhai  et  al.  calculated  the  activation  energies  for  vacancy 
migration  along  different  paths  in  praseodymium  doped  ceria  [38] 
and  gadolinium  doped  ceria  [36].  The  former  can  be  extended  to 
represent  the  cases  where  the  ionic  radius  dopant  ion  is  larger  than 
that  of  the  host  ion;  while  the  latter  can  be  extended  to  repre¬ 
sent  the  case  where  the  ionic  radii  of  the  dopant  and  the  host  have 
similar  values.  According  to  their  calculation  results,  in  both  cases 
the  preferred  migration  pathway  for  oxygen  vacancy  diffusion  is 
between  1  NN  and  2NN  sites.  For  example,  in  yttria  doped  ceria,  the 
oxygen  vacancy  is  located  at  the  1  NN  position.  For  migration,  it  will 
jump  to  2NN  position  first,  and  then  it  is  of  high  possibility  that  it 
jumps  back  to  the  INN  site  because  the  jump  from  2NN  to  INN 
has  lower  activation  energy  compared  with  the  jumps  from  2NN  to 
2NN  or  from  2NN  to  3NN  sites.  The  back-and-forth  jump  between 
two  sites  contributes  to  the  localized  conduction. 

In  summary,  the  localized  relaxation  is  from  the  orientation  of 
the  defect  associates  under  the  electric  field  and  the  back-and-forth 
jump  between  two  sites  of  a  free  oxygen  vacancy. 

4.3.  Composition  dependence  of  the  ionic  conductivity  and  the 
activation  energy 

In  Fig.  6,  a  “V-shape”  variation  of  the  bulk  conductivity  as 
a  function  of  x  is  observed.  Similar  trend  of  the  composition 


dependence  of  the  electrical  conductivity  in  Zr02-Ce02-Y203  solid 
solutions  has  been  reported  by  Ananthapadmanabhan  et  al.  [24]. 
They  attributed  this  V-shape  variation  to  the  increased  scatter¬ 
ing  of  the  oxygen  ions  by  the  substitution  of  Zr  by  Ce  atoms  or 
Ce  by  Zr  atoms,  which  is  analogical  to  the  mechanism  of  enhanced 
electrical  resistivity  in  an  alloy  system.  This  interpretation  is  uncon¬ 
vincing  to  some  point  that  the  hopping  mechanism  in  an  ionic 
conductor  is  totally  different  from  the  scattering  mechanism  of 
electrons  in  metals  and  alloys.  Kawamura  et  al.  [20]  also  found 
the  “V-shape”  composition  dependence  of  the  conductivity  in  CaO 
doped  Zr02-Ce02  solid  solutions,  however  they  did  not  give  any 
further  explanations  on  the  results.  In  the  following  paragraphs 
the  composition  dependence  of  the  conductivity  will  be  discussed. 

Within  the  studied  temperature  range,  the  electrical  conduction 
in  the  solid  solutions  is  purely  ionic.  Therefore,  the  variation  of  the 
bulk  conductivity  as  a  function  of  the  composition  is  either  caused 
by  the  change  of  the  oxygen  vacancy  concentration  or  its  mobility. 
Here  we  divide  the  “V-shape”  variation  as  a  function  of  the  com¬ 
position  in  Fig.  6  into  two  regions:  the  Ce02-rich  (x>0.5)  and  the 
Zr02-rich  (x<  0.5)  regions  and  discuss  them  separately. 

The  increase  of  the  ionic  conductivity  with  the  increase  of  x  in 
the  Ce02-rich  region  is  easier  to  understand  to  some  extent.  In  the 
fluorite  structure,  the  ionic  conduction  is  caused  by  the  migration 
of  oxygen  vacancies  through  channels  formed  by  the  neighbour¬ 
ing  cations.  As  discussed  by  Ananthapadmanabhan  et  al.  [24]  and 
Tsoga  et  al.  [20],  the  radius  of  the  free  channel  increases  with  the 
lattice  parameter.  Larger  channel  radius  is  favourable  for  the  charge 
carriers  to  migrate.  Since  the  amount  of  Y203  is  set  constant  for  all 
the  compositions,  the  concentration  of  oxygen  vacancies  remains 
almost  the  same  (the  slight  increase  in  the  lattice  volume  leads 
to  a  slight  decrease  of  the  defect  concentration,  but  this  effect  is 
negligible  [24]).  Therefore,  the  increase  of  ionic  conductivity  in 
the  Ce02-rich  region  with  increasing  x  can  be  attributed  to  the 
increased  oxygen  vacancy  mobility  through  a  larger  channel. 

However,  in  the  Zr02-rich  region  (x<0.5),  the  increased  chan¬ 
nel  radius  fails  to  explain  the  decreased  ionic  conductivity  with 
increasing  x.  The  possible  reason  for  the  decreased  conductivity  in 
the  Zr02-rich  region  with  increasing  Ce02  mole  ratio  is  the  for¬ 
mation  of  complex  defect  associates  or  clusters.  It  is  well  known 
that  8mol%  YSZ  has  the  highest  electrical  conductivity.  Further 
increases  in  dopant  concentration  decrease  the  conductivity  in  YSZ 
because  of  the  formation  of  defect  associates,  which  bind  the  oxy¬ 
gen  vacancies  to  the  yttrium  ions,  making  the  oxygen  vacancies 
unavailable  for  conduction  [25].  The  introduction  of  Ce02  into  YSZ 
has  similar  effect:  when  Ce02  is  introduced  into  Zr02,  additional 
oxygen  vacancies  can  be  generated  due  to  the  significant  differ¬ 
ence  of  the  ionic  radius  between  Ce4+  and  Zr4+  (which  is  called 
“size  effect”)  by  the  following  reaction  [39]: 

Ce02^Ce*r  +  2Vq  +  20". 

The  oxygen  vacancies  generated  by  Ce02  addition  makes  the  for¬ 
mation  of  (Y^Vq  )  becomes  more  significant,  and  possibly  leads  to 
more  complex  defect  associates  (Y^Vq  ,  which  has  been  found 

in  9  mol%  Y203  doped  Zr02  [40].  The  formation  of  defect  associates 
has  stronger  effect  over  the  increase  in  the  lattice  parameter,  and 
consequently  leads  to  a  decreasing  conductivity. 

To  sum  up,  the  “V-shape”  variation  of  the  bulk  conductivity  as 
a  function  of  x  in  the  solid  solutions  is  a  competition  between 
two  factors.  The  first  one  is  the  formation  of  defect  associates 
which  limits  the  mobility  of  the  oxygen  vacancies  and  therefore 
decreases  the  conductivity;  the  second  one  is  the  increase  of  lat¬ 
tice  parameter  which  enlarges  the  free  channel  for  oxygen  vacancy 
migration  and  therefore  increases  the  conductivity.  In  the  Zr02-rich 
region  (x  <  0.5),  the  formation  of  defect  associates  overshadows  the 
effect  of  the  increased  lattice  parameter,  therefore  the  conductiv- 
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Fig.  9.  (a)  Schematic  of  the  defect  associates  profile  (black  dash  line)  and  lattice  parameter  (red  dot  line)  as  a  function  of  the  composition  (x).  The  black  solid  line  illustrates 
the  conductivity  decrease  caused  by  increasing  defect  associates;  the  red  solid  line  illustrates  the  conductivity  increase  due  to  larger  lattice  parameter;  the  blue  solid  line 
shows  the  competitive  effect  of  the  two  factors  on  the  bulk  conductivity  of  the  solid  solutions,  (b)  Schematic  of  the  competitive  effect  of  lattice  parameter  and  the  defect 
associates  on  the  activation  energy  for  the  bulk  conduction.  (For  interpretation  of  the  references  to  color  in  this  figure  caption,  the  reader  is  referred  to  the  web  version  of 
the  article.) 


ity  decreases  with  increasing  x;  in  the  Ce02-rich  region  (x  >  0.5),  the 
increase  of  lattice  parameter  is  the  dominant  factor  which  leads  to 
an  up-going  trend  on  the  conductivity.  A  schematic  of  the  defect 
associates  profile  and  lattice  parameter  as  a  function  of  x  and  their 
competitive  effect  on  the  bulk  conductivity  of  the  solid  solutions  is 
shown  in  Fig.  9(a). 

The  defect  associates  and  the  lattice  parameter  also  influence 
the  activation  energy  for  the  bulk  conduction,  as  illustrated  in 
Fig.  9(b).  It  is  accepted  that  the  activation  energy  for  the  oxygen 
diffusion  (oxygen  vacancy  diffusion  in  an  opposite  direction)  in  the 
intermediate  temperature  range  is  comprised  of  migration  energy 
for  a  free  oxygen  vacancy  (EM)  and  the  dissociation  energy  (E0) 
to  extricate  an  oxygen  vacancy  from  the  defect  associates.  As  dis¬ 
cussed  above,  the  increase  of  lattice  parameter  leads  to  more  free 
space  and  makes  it  much  easier  for  an  oxygen  vacancy  to  migrate, 
and  therefore  a  lower  EM.  Eo  is  dependent  on  the  concentration  of 
defect  associates.  More  defect  associates,  more  energy  required  to 
dissociate  them.  Therefore  E0  varies  as  a  function  of  x  in  the  same 
way  as  the  concentration  of  defect  associates.  Since  EM  has  higher 
values  than  E0  [41  ],  in  the  Zr02-rich  region,  the  increased  E0  slows 
down  the  decrease  of  EM,  therefore  the  activation  energy  shows 
quite  smooth  variation;  while  a  rapid  drop  of  the  activation  energy 
can  be  observed  in  the  Ce02-rich  region  because  EM  and  E0  both 
decrease  as  a  function  of  x.  The  experimental  observation  in  Fig.  6 
is  in  agreement  with  the  above  analysis. 

From  the  above  analysis  it  can  be  concluded  that  the  “V-shape” 
variation  of  the  bulk  conductivity  with  increasing  Ce02  mole  ratio 
is  because  the  ionic  radius  of  Ce4+  is  significantly  larger  than  Zr4+. 
Therefore,  the  ionic  radius  of  the  tetravalent  dopant  determines 
the  composition  dependence  of  the  ionic  conductivity  of  the  solid 
solutions.  It  can  be  predicted  if  Zr  is  substituted  by  another  tetrava¬ 
lent  element  with  the  similar  ionic  radius  with  Zr4+  (e.g.,  Hf), 
“V-shape”  composition  dependence  of  the  bulk  conductivity  may 
not  be  observed. 

In  order  to  testify  the  above  prediction,  another  series  of  binary 
solid  solution  [(ZrO2)i_x(HfO2)x]0.87(Y2O3)0.i3  (0<x<0.8)  has 
been  investigated.  The  bulk  conductivity  and  the  activation  energy 
were  obtained  using  the  same  method  as  the  Zr02-Ce02-Y203 
series.  The  ionic  radius  of  Hf4"*"  is  0.083  nm,  almost  the  same  as  Zr4+ 
(0.084  nm),  thus  the  substitution  of  Zr4+  by  Hf4"1"  will  not  cause  the 
re-arrangement  of  oxygen  vacancies  (no  “size  effect”).  Besides  at 
given  yttria  concentration,  yttria-stabilized  hafnia  has  lower  ionic 
conductivity  and  higher  activation  energy  than  YSZ  [42],  therefore 
we  expect  a  monotonously  decrease  of  the  bulk  conductivity,  along 


with  a  monotonously  increase  of  the  activation  energy  as  a  func¬ 
tion  of  x.  As  shown  in  Fig.  10(a)  and  (b),  both  the  variations  of  the 
bulk  conductivity  and  the  activation  energy  as  a  function  of  x  are 
in  accordance  with  the  expectations.  No  “V-shape”  composition 
dependence  of  the  bulk  conductivity  is  observed,  which  confirms 
the  prediction  and  from  the  other  side  supports  the  validity  of  the 


Fig.  10.  Composition  dependence  of  (a)  bulk  conductivity  and  (b)  activation  energy 
of  the  [(Zr02  )i_x(Hf02  )x]o.87(Y203  )0.i3  (0  <x  <  0.8)  solid  solutions  at  various  temper¬ 
atures. 
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above  explanation  of  the  “V-shape”  composition  dependence  of  the 
bulk  conductivity  in  Zr02-Ce02-Y203  solid  solutions. 

5.  Conclusions 

In  this  paper,  the  bulk  conduction  and  relaxation  of  the 
[(Zr02)i_x(Ce02)x]0.92(Y203)o.o8  (0<x<l)  solid  solutions  were 
studied  using  impedance  spectroscopy  at  intermediate  tempera¬ 
tures  (200-500  °C).  The  main  conclusions  are:  (1)  the  impedance 
spectra  of  the  solid  solutions  of  all  the  compositions  show  similar 
characteristic  of  the  well  established  impedance  spectra  of  YSZ;  (2) 
the  relaxation  frequency  (fz»  and  fM" )  is  proportional  to  the  bulk 
conductivity.  Comparing  the  relaxation  frequency  can  be  a  simple 
method  to  rank  the  bulk  conductivity  of  the  solid  solutions;  (3)  the 
electrical  response  of  the  solid  solutions  for  all  the  compositions 
under  the  ac  electrical  field  shows  both  long-range  and  localized 
relaxations.  The  localized  relaxations  are  from  the  orientation  of 
the  defect  associates  under  the  electric  field,  as  well  as  the  back- 
and-forth  jumping  of  a  free  oxygen  vacancy  between  the  two  sites 
through  which  the  diffusion  path  has  the  lowest  activation  energy; 
(4)  the  bulk  conductivity  has  a  “V-shape”  variation  as  a  function  of 
the  composition:  it  decreases  with  an  increase  of  the  mole  ratio  of 
Ce02  (x),  reaching  the  lowest  value  when  x  =  0.5  (Zr/Ce  =  1 ),  and  then 
goes  up  with  further  increase  ofx.  The  decrease  of  the  bulk  conduc¬ 
tivity  with  increasing  x  in  the  Zr02-rich  region  (x  <  0.5)  is  due  to  the 
formation  of  the  defect  associates  which  limits  the  mobility  of  the 
oxygen  vacancies,  while  the  increase  of  the  bulk  conductivity  with 
increasing  x  in  the  Ce02-rich  region  (x>0.5)  is  due  to  the  increase 
of  the  lattice  parameter  which  enlarges  the  free  channel  for  oxygen 
vacancy  migration.  These  two  factors  also  influence  the  activation 
energy  for  the  bulk  conduction;  (5)  the  ionic  radius  of  the  tetrava- 
lent  dopant  determines  the  composition  dependence  of  the  ionic 
conductivity  of  the  solid  solutions.  When  doping  YSZ  with  other 
tetravalent  dopant  with  similar  ionic  radius  with  Zr4+,  e.g.,  Hf4+, 
such  “V-shape”  composition  dependence  of  the  bulk  conductivity 
cannot  be  observed. 
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